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Abstract 

A systematic procedure to consistently formulate a field theoretical, QCD 
bound state problem with a fixed number of constituents is outlined. The 
approach entails applying the Hamiltonian flow equations, which are a set of 
continuous unitary transformations, to a QCD motivated Hamiltonian with 
a confining interaction. The method is developed in detail for gluodynam- 
ics in the Coulomb gauge to obtain an effective block-diagonal Hamiltonian 
appropriate to a reduced Fock space with fixed number of dynamical glu- 
ons. Standard many-body techniques are used to numerically diagonalize this 
Hamiltonian in a constituent two gluon Fock space. The calculated gluon 
condensates and glueball masses are in good agreement with QCD sum rule 
and lattice results. 

PAC number(s): 12.39.Mk, 12.39.Ki, ll.10.Ef, ll.10.Gh 
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I. INTRODUCTION 



Quantum Chromodynamics (QCD) is generally regarded as the theory governing 
hadronic structure and interactions. This broad acceptance follows, in part, from the suc- 
cess of high energy perturbative calculations where asymptotic freedom prevails. However, a 
significant gap still exists in implementing QCD in the non-perturbative, low energy regime 
where physical observables have only been described by lattice calculations, QCD sum rules 
and phenomenological models. Our study addresses this shortcomming and details a new 
approach which directly incorporates canonical QCD at high energies and phenomenology, 
including confinement, ^from low energy, to produce a realistic effective Hamiltonian appro- 
priate for bound state investigations. 

Because the relevant hadronic degrees of freedom do not clearly follow from the canonical 
QCD Hamiltonian, we look to phenomenological models, which are reasonably developed, for 
initial guidance. Of particular interest is gluonic structure, the focus of this work, which has 
been described by vastly alternative approaches: strings [|IJ, flux tubes 0, stochastic gluon 
configurations [|| and constituents [f| confined by a bag || or potential ||[7|]. In this paper we 
continue to build from our previous work |7]-|9| within the constituent picture, again viewing 
both quarks and gluons as dressed partons subject to an effective interaction. We assume 
this interaction is dominated by a confining potential whose explicit form is not essential, 
provided it is long ranged and generates bound states. Our effective Hamiltonian also 
contains perturbative interactions, including radiative corrections, governing the ultraviolet 
(UV) region. These terms are directly and unambiguously obtained from the leading order 
coupling constant expansion of the canonical Hamiltonian. 

With these preliminaries we now address the essence of this paper, a systematic method 
for consistently solving a generic, field theoretical QCD motivated Hamiltonian. The method 



we present is based upon the flow equations suggested by Wegner |1(J and applied to QED 
within the light-cone framework [ 0], l2| . Because a field theoretical solution involves diag- 
onalization in an infinite Fock space with arbitrary number of quanta, an exact solution 
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is not possible. Instead we seek a finite Fock space formulation and begin by dividing the 
complete Fock space into two components, a tractable P subspace spanned by states with a 
small number of quanta and the remainder Q = 1 — P. This leads to a Hamiltonian matrix 
to be diagonalized of the form 



H 



( PHP PHQ ^ 



V 



QHP QHQ 

Next we use the flow equation method to block diagonalize this matrix yielding the effective 
Hamiltonian 



H. 



eff 



V 



X2) 



J 



PHeffP 

QH eff Q 

One can then separately diagonalize the two blocks which are now uncoupled. It is important 
to note that the flow equation approach is significantly different than the Tamm Dancoff 
approximation (TDA) even though they both may use the same P space. In the TDA the 
off-diagonal blocks PHQ and QHP are eliminated simply by truncation, whereas in the 
flow equation method they are eliminated by transformation to PH e ffP and QH e f/Q. 
In the flow equation approach the Hamiltonian is unitary transformed continuously 



H(l) = U{l)H(0)U\l) 



;i.3) 



U(l)=T l exp ^J\(l')dl^j , (1.4) 
where / is the continuous (flow) parameter (1 = corresponds to the initial Hamiltonian, 



Eq. (1.1)), Ti the / ordering operator and rj the generator of the transformation. Wegner's 
flow equations are then an expression of this unitary transformation in differential form 

dH(l) 



dl 



[ V (l),H(l)] 



;i.s) 



with the specific generator choice 



17(0 = ^(0^(0] 



;i.6) 



where Hd is the diagonal part of H and H e ff = Hd(l — > oo). Using this generator, it is 
shown in section II how the off-diagonal terms QHP, which connect Fock states with different 
particle number, can be eliminated as I — > oo and included in the block diagonal effective 
interaction. Hence coupling to the Q space is incorporated into the P space and represented 
by new terms in the effective Hamiltonian. Once PH e ffP is determined, the bound state 
problem can be investigated in detail by direct diagonalization. In general, PH e ffP will 
become more complicated as the P space becomes smaller which may necessitate additional 
approximations in obtaining the effective Hamitlonian. Hence the strategy it to find a 
balance between size of the P space, which is determined by computational effort and 
resources, and level of approximations. In the current work we will apply this approach to 
gluonia and restrict the P space to Fock states with at most two dynamical (dressed) gluons. 

In the next section we review the flow equation approach for a general Hamiltonian. It 
is important to note, however, that the solution of these equations is quite involved. This 
is because the equations are nonlinear and coupled and also because the choice of generator 
needed to decouple the P and Q spaces generates complicated many-body interactions. For 
most realistic applications, therefore, the flow equations will be implemented iteratively. 
This has been demonstrated in a variety of investigations such as condense matter problems 
|I3| , |14|| , electron-phonon Jl5| and positronium treatments, the latter two involving 



perturbative expansions in coupling constant. Here we adopt a slightly different perturba- 
tive scheme and invoke the weak coupling assumption commonly used in nuclear structure 
many-body calculations. There one introduces a mean field, or Hartree-Fock component 
of the interaction, and expresses the Hamiltonian interaction as a sum of this term and a 
residual force, the difference between the original potential energy and the mean field. An 
approximate effective interaction is then generated by treating the residual interaction per- 
turbatively. The physical nuclear many-body states are finally obtained by diagonalizing this 
effective interaction in the P space spanned by eigenf unctions of the mean field Hamiltonian. 
The success of the method is determined by choice of model space P and phenomenological 
interactions. Although our problem is not a central field situation, we nevertheless can uti- 



lize this concept since a dominant mean field is present in the form of the phenomenological 
confining interaction. This directly governs the quasi-particle (dressed) gluon basis functions 
which in turn form the many-body basis spanning the P space. We identify our perturbative 
interaction as the difference between the original QCD interaction and the phenomenolog- 
ical confining potential and approximate this residual interaction at high energies by the 
leading order canonical Hamiltonian interaction. Again the overall utility of this method 
will depend significantly upon choice of P space and phenomenological input. For example 
in our pure gluon application we assume that three and higher gluon Fock states only couple 
weakly to the P space. Because the dressed gluon mass is of order 1 GeV, which sets the 
scale for a large energy denominator in a perturbative expansion, this approximation should 
be reasonable. However, it may develop in quark sector applications that the weak coupling 
approximation is not sufficient for low Fock space components such as exotic 4 quark (qqqq) 
and hybrid (qqg) states. In such cases the P space must be expanded to explicitly include 
these higher components and this is a subject for future investigations. 

Finally we mention the issue of renormalization which can also be accomodated by the 
flow equations. Because the complete Fock state expansion contains states of arbitrarily 
large energies, UV divergences occur for certain matrix elements. However, since QCD is a 
renormalizable theory, these UV diverging terms can be isolated, regularized and then by 
introducing counterterms, renormalization can be achieved perturbatively order by order in 
the coupling constant Pjl6|j. As demonstrated in the following sections, the flow parameter, 



actually Z -1 / 2 , plays an additional role as an UV cut-off in the flow equation elimination 
of the off-diagonal terms. Hence, the flow equations embody the renormalization concept 
of Wilson wherein the high energy modes are integrated out in the path integral represen- 
tation yielding a low energy effective action. Our method is also analogous to the Glazek 
and Wilson similarity transformation scheme [|nj with a continuous cut-off function ||. In 
particular, as shown in section IV. A, we obtain a well behaved gap equation with canonical 
counterterms cancelling the leading UV divergence. This is one reason why we do not incor- 
porate the phenomenological confining terms in the implementation of the flow equations 



since they are not renormalizable. Further and related, including the confining potential in 
the flow generator would produce additional effective interaction terms that are still phe- 
nomenolgical. Since an initial phenomenological interaction is already present, the model 
framework remains unchanged and we therefore bypass this step. 

This paper is divided into five sections and three appendices. In the next section we 
summarize the general flow equation method, highlighting the key issues. In section III we 
apply the flow equations and specify the explicit form of the model QCD Hamiltonian. In 
addition to the weak coupling motivation discussed above, we also wish to disentangle all 
complexities from chiral symmetry breaking and thus focus on pure gluodynamics. This 
section also developes and presents our final effective Hamiltonian. Numerical results are 
then given in section IV, with the gap equation, gluon constituent mass and condensate 
analyzed in section IV.A and glueball bound states and masses in section IV.B. Finally, 
results and future work are discussed in the concluding section along with an extensive 
appendix containing full details of our comprehensive treatment. 



II. THE FLOW EQUATION METHOD 



We now summarize the important elements of flow equations for a general Hamiltonian. 



We follow Wegner |T0[ who first formulated this approach to construct an effective bound 
state Hamiltonian appropriate for field theories. Since this method does not depend on the 
specific system or detailed nature of the Hamiltonian, it is beneficial to first conceptually 
outline the general idea before our application in section III to a complicated field theory. 



See ref. [0 for complete details. 

Continuing with Eqs. (1.3-1.6), the first step is to divide the Hamiltonian, H(l), into 
diagonal and "rest", or off-diagonal, components, H d (l) and H r (l), respectively 



H d (l) 



( PH{l)P ^ 



V 



QH{l)Q 



(2.1) 



6 



H r (l) = H(l) - H d (l) = 



( PH(l)Q ^ 



(2.2) 



\QH(l)P J 

It is precisely the off-diagonal part, which connects Fock states with different particle num- 
ber, that we wish to remove. Next, and most important, is the choice of generator. As 
shown by Wegner our block diagonalization goal can be acheived by choosing 

V (l) = [H d (l),H(l)\, (2.3) 



which is off-diagonal by direct evaluation 



( P V {l)Q ^ 



(2.4) 



K Q V (l)P } 

with Pr/Q = PHPHQ — PHQHQ and similar result for Qr/P upon interchange of P and 
Q. The flow equation, Eq. (1.5), can then be written for the important sector Hamiltonians 
PHP and PHQ 

%PHP = PriQHP - PHQr]P , (2.5) 
dl 

jfHQ = PrjQHQ ~ PHPrjQ . (2.6) 

The essence of the flow equation approach is that with Wegner's generator choice the absolute 
value of the off-diagonal terms, PH(l)Q, QH(l)P, monotonically decreases to zero as the 
flow parameter tends to infinity. This can be directly shown by taking the trace of their 
product and computing the change with I yielding, using EQ. (2.6), a negative derivative 

^-TtPHQHP = Tr(Pr]Q(QHQHP - QHPHP)) 
+ Tr((PHPHQ - PHQHQ)Qr 1 P) 
= 2Tr (PrjQrjP) 

= -2Tr(P??Q(P?7Q)t) < 0, (2.7) 

where the last step follows from the antihermiticity of PrjQ. Since the trace of PHQ HP 
is positive definite, Prj(l)Q — > in the limit I — > oo. Hence, H d commutes with H and the 



effective Hamiltonian, H e ff, approaches complete block diagonal form with increasing flow 
parameter. 

We can formally solve the flow equations to provide further detailed insight and also 
the starting point for our gluodynamics application in the next section. The Hamiltonian, 
H = H-f ree + H 1 , consists of free (non-interacting), H free , and interacting, H 1 , components. 
We then work in a fixed basis spanned by eigenfunctions of the free sector Hamiltonians 
pjjfreep an( j QjjfreeQ w ^ eigenvalues E p and E q , respectively. Here the indices p and q 
denote all states (both Fock and energy) in the respective P and Q spaces. Representing 
the matrix elements of transformed Hamiltonians PH(l)P, PH(l)Q, etc. in this basis by 
h pp '(l), h pq {l), ... the flow equations take the form 

^f" = 12 (jlpqhqp' - KqVqp') , (2-8) 

^ = 12 Vpq'hq'Q ~ 12 Kp'Vp'q , (2-9) 
q' p' 

Vpq — 12 hpp'K'q ~ 12 Kq'hq'q > (2-10) 
p' q' 

where all matrix elements are functions of I. The initial conditions (flow parameter / = 0) for 
the matrix elements of the sector Hamiltonian PH(0)P are, h pp '(0) = Ep(0)5 pp '+5E p T (0)5 pp i + 
hp p ,(l — S P pf) = Ep(0)5 pp i + hp^lp^p' , and for the Hamiltonian PH(0)Q, h pq (0). Here 5E^ and 
hppi\p^p' are the diagonal and off-diagonal matrix elements of the interaction Hamiltonian 
H 1 in the P space. We now separate the diagonal and off-diagonal matrix elements in Eqs. 
(2.8-2.10) 

dE^ 

~ = 12 (VpqKp - h pqVqp) , (2-H) 



dl 

dhL 



PP 



dl 



12 (VpqKp' - KqVqp') , (2-12) 



= " {K - Ei) V p q + £ Vpq'K,, ~ 12 h^Vp'q , (2.13) 

q'+q p'+p 

Vpq = (#J - E T q ) h pq + ]T hI PP ' h p'q - 12 Kq'K'q ■ ( 2 - 14 ) 

p'+p q'^q 

This system of coupled equations can be approximately decoupled by assuming that the off- 
diagonal PHQ sector matrix elements are small as compared to diagonal PHP elements, 



i.e. \h pq (0)\ <C \h pp /(0)\ (from the flow equations this is also valid for any finite /). This 
approximation is consistent with our weak coupling concept discussed earlier and is valid for 
our gluonic application which entails a nonperturbative confining potential in the diagonal 
sector Hamiltonian and a weak, perturbative residual interaction in the off-diagonal. 

We now solve the flow equations perturbatively, iterating h pq . We do not treat h pp > 
perturbatively, but note it depends very weakly upon I, being suppressed to second order in 
h pq since 'Jiom. Eqs. (2.11,2.12), dE p /dl ~ 0(h pq ) and dh pp ,/dl ~ 0(h pq ). Hence, our leading 
iteration in Eqs. (2.13,2.14) uses E* p (l) ~ ££(0) = E p (0) + SE^O) and h^(l) ~ h^(0). We 
could systematically iterate Eqs. (2.11-2.14) carrying all terms, however, it is simpler and 
more direct to further approximate Eqs. (2.13,2.14) by recognizing that the two off-diagonal 
sums over p' and q' will tend to be suppressed by cancellation. To appreciate this point, 
consider the leading order iteration of Eq. (2.14) in Eq. (2.13) which generates terms of 
the type, J2 P > hp\,,(0)ri P > q (l) ~ J2 p >(E p ,(0) - E q (0))h ppl (0)h p / q (l), with a similar expression for 
the Q-space sum, J2 q > Vpq'i^hLJO)- Here the sum is over Fock state components and all 
momentum (energy) for each Fock state. Next we make the key observation that although 
the Q space has in general an infinite number of Fock states, the Fock sum over q' is restricted 
to at most two adjacent sectors since h qq > is a two-body interaction which can change particle 
number at most by two. Hence the P and Q space "volumes" in particle number (and energy 
for a given particle number sector) are the same. The two sums are therefore equivalent to 
the product of an averaged matrix element for each space with a common volume for both 
spaces. Because the matrix elements have random phases, in the limit of large numbers 
of energy states, the two sums will tend to cancel and thus be suppressed relative to the 
leading diagonal terms (p — p',q — q'). 

Therefore in leading order the generator is 

Tfo® = ~ E}(0) - E' q (0) (2 ' 15) 
and the effective Hamiltonian PH e ffP is given by h p J = h pp >(oo) 
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roc ( dh vS) h ,m u (n dh «p'( l ) \ 

W - KW* + *^(o)U - i - 4^s) ' (2 - 16) 

having off-diagonal elements 

MO = V(0)exp (-/(^(0) - ^(O)) 2 ) . (2.17) 



which can be used for the next iteration. The above process can then be repeated, replacing 
Ep(0), h pp '(0) on the right hand side of the flow equations with the diagonal and off-diagonal 
matix elements of h e JJ , respectively, and continued until convergence. Diagonalizing h e p J 
then produces the final eigenvalues E%H. We now apply these results to gluodynamics. 



III. EFFECTIVE HAMILTONIAN FOR GLUODYNAMICS 

In constructing our QCD motivated Hamiltonian, Hqcd-, we have incorporated phe- 
nomenological elements of previous, related constituent studies |7|||. Starting with the 
exact Coulomb gauge QCD Hamiltonian (pure gluodynamics) fll?! , Hqcd = Hf ree + Hq CD , 
containing free and QCD interacting terms, we introduce a confining phenomenological in- 
teraction, H p hen, to permit isolation of the nonperturbative and perturbative Hamiltonians, 
Hjsrp and Hpx, respectively 

Hqcd = H NP + H PT = H NP + B.]. an = Hq C d > (3.1) 
H NP = Hf ree + Hphen , (3.2) 

H PT = Hq CD - H phen = Hj, an . (3.3) 

The soluable, phenomenological interaction H p h e n is introduced to provide a weaker, residual 
interaction Hpt at all energy scales. Following quark models, we assume that, in the con- 
stituent (quasiparticle) basis, H NP predominantly governs the bound state bulk properties 
and should be treated nonperturbatively. We further assume that the residual interac- 
tion can be approximated at high energies by the canonical QCD interaction Hamiltonian 
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Hpt = H^ an (Auv oo) and is amenable to perturbation theory. Note that the residual in- 
teraction describes quantum fluctuations and may contribute UV divergent terms requiring 
regularization with cut-off parameter Auv- 

The goal is to use the flow equations to scale the residual interaction to low energies 
to generate an effective interaction, which when added to H^p, can be diagonalized non- 
perturbatively for bound states. In this process renormalization will also be achieved by 
introducing second order 0(g 2 ) canonical counterterms which remove the leading UV diver- 
gences. With an assumed weak residual interaction, the entire procedure can be conducted 
perturbatively. There is no over representation or double counting, since, as shown below, 
the phenomenological interactions determine the infrared (IR) behavior while the residual 
interactions, which are free from IR divergences, describe the UV part. 

We now further detail the Hamiltonian. In the Coulomb gauge the physical degrees of 
freedom are the transverse gauge fields A a T a (with color matrices T a and sum over index 
a) and their conjugate transverse momenta II. The free Hamiltonian is 

H free = Tr J rfx (n 2 (x) + B\ (x)) , (3.4) 

where = V x A is the abelian component of the color magnetic field. The nonabelian 
contribution, B 2 (x) — B^(x), is included in Hpt (see below). Here B = B a T a is the com- 
plete magnetic field having components B% = e^V ' jA% + ^eijkf abc ' A h -A c k , with bare coupling 
constant g. As in the phenomenological Hamiltonian is taken to be 

H phen = ~J dxcfyp a (x)V L (|x - y|)p a (y) , (3.5) 

with color charge density (only gluon component) p a (x) = / abc A 6 (x)II c (x) and linear con- 
fining potential 

\4(|x-y|) = a|x-y|, (3.6) 

with string tension, a = 0.18 GeV 2 , in accordance with lattice calculations and Regge phe- 
nomenology. The approximate residual interaction is given by the sum, Hpt = H^a + He, 
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with the nonabelian part including three-, H^ g , and four-gluon, H^g, interactions (see Ap- 
pendix A) 

H NA = Tr J rfx (B 2 (x) - Bi(x)) = H 3g + H 4g , (3.7) 
and the Coulomb interaction 

Hc = -\J dxdyp a (x)V c (|x - y\)p a {y) , (3.8) 

lfc(|x-y|) = -C^^, (3.9) 

having a s = ^- and adjoint Casimir operator C a dj = (N% — l)/2iV c (here N c = 3 colors). An 
advantage of the Coulomb gauge is the straightforward implemention of phenomenological 
potentials. We define the instantaneous interaction, H^ + c, by the combination 

Hl+c = H phen + He , (3.10) 

containing the sum, Vl + Vc, of linear and Coulomb potentials. 

Because of the phenomenological interaction, the trivial, perturbative vacuum is not the 
minimum ground state for the QCD motivated Hamiltonian, Hqcd- Instead, we introduce 
a trial nonperturbative vacuum state |0) which can be determined variationally, analogous 
to BSC type vacuum |7|]8|,|l8-|2~0| studies. The Fock space is constructed ^from this vacuum 
using quasiparticle operators a"(k), a%' (— k) which appear in the field operator expansions 

• 4 ° (x) = /(^7fc K(k) + af( - k)|e ' kx 

n?(x) = -iJ ^^[<(k) -<* (-k)]e ta , (3.11) 

with a*(k)|0) = 0. Here is the variational gluon energy from the gap equation (see below), 
obtained by minimizing the ground state (vacuum) energy. The canonical commutation 
relation is 

[at(k),af(k')} = (27r)W 3 )(k-k')A,(k), (3.12) 
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where the gluon operators a"(k) = J2x=i,2 e «(k, A)a a (k, A) are transverse, i.e. k- a a (k) = 
k • a a ^(k) = 0, and -Dy(k) is a polarization sum 

D i:j (k) = e *(X A)ej(k, A) = - kikj , (3.13) 

A=l,2 

with unit vector component fcj = fcj/A; and fc, • Dy(k) = 0. 

The complete QCD motivated Hamiltonian, Eq. ( |3.1| ), can now be expressed in second 
quantized form and normal ordered with respect to the trial vacuum. In addition to two- 
body interactions, normal ordering leads to condensate (vacuum expectations) and one-body 
operators, both having contributions from the instantaneous and four-gluon terms. The final 
expressions are complicated and are summarized in Appendix A. 

Finally the flow equations are applied to Hqcd = Hf ree + Hl + c + H 3g + H^. Seperating 
the diagonal, Hj, particle number conserving and off-diagonal, H r , particle number changing 
parts yields 

H d =(Hfl) +«) , 

\ J / ran V ' J nn V a / nn 

^=«) +« e ) +(4 2) ) +«) , (3-14) 

V a / ram \ J / ram V / nm V M / nm 

where nn (nm) indicates particle number conserving (changing) matrix elements. Here and 
throughout, the superscript denotes the order (power) in the bare coupling constant. Note 
that the triple gluon Hamiltonian, H^, which is first order, only contributes off-diagonally, 

(2) 

while the four-gluon Hamiltonian, H\ g , which is second order, contributes to both and 
H r . Also, as discussed previously, the linear confining potential, which is zeroth order, is 
restricted to diagonal contributions. All terms are detailed in Appendix A along with the 
free Hamiltonian operator, containing the kinetic (K,K^°>) and condensate (0,0^) terms, 
which contributes to the diagonal part in lowest order, HfJ ee = + 0^°\ and to the 
off-diagonal part in second order, nf} ee = (K — K^) + (O — O^). 

Having identified the diagonal and rest parts we now construct the generator, rj, and 
solve the flow equations for the first two leading orders. The leading order 0(g) (in coupling 
constant) flow equation is 
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V ^(l) = [K(°\l),H^(l)]. (3.15) 
Here the kinetic and triple-gluon vertex terms are given (see Appendix A) respectively by 

^ (0) (0 = / (|j^( k ; 0«f(k)<(k) , (3.16) 
and in symmetrized form 

«u = ^/-/ (n^) 



x: 



p (ki, k 2 , k 3 ; /)<(k 1 )a;(k 2 )a c fc (k 3 ) + 3^(k x , k 2 , k 3 ; Z)a? t (-k 1 )a5(k 2 )o^(k 3 ) + h.c. : , (3.17) 



with shorthand notation lo\ = cj(ki), etc. and 1^ defined by 



r yfc (ki, k 2 , k 3 ) = i [(fei - k 3 )j5 ik + (fc 2 - + (fc 3 - fca)^] , 



satisfying 



rijk(-ki, -k 2 , -kg) = -T ijk (k 1 , k 2 , kg) , 

T jik (k 2 , ki, kg) = -r iifc (ki, k 2 , k 3 ) . (3.18) 

Notice that in implementing the flow equations new, effective coupling constants are gener- 
ated (see below), <7o(ki, k 2 , k 3 ; I), (?i(k 1; k 2 , k 3 ; I), which like the energy u;(k; /), depend upon 
the flow parameter I. The coupling constants also aquire an effective energy /momentum 
dependence corresponding to a given Fock sector (different Fock sector operators flow differ- 
ently in energy, indicated by and 1 subscripts here and below). In leading order we omit 
the dependence of gluon energy on I in the triple vertex. The leading order generator, Eq. 
(E3), is 



^(0 = ^f abc I (ft Pk] (2vr) 3 ^(Ek m )-^=r^( kl ,k 2 ,k 3 ) (3.19) 



2^/2 

x: 



77o(ki, k 2 , k 3 ; /)<(k 1 )a J b (k 2 )4(k 3 ) + 3^, k 2 , k 3 ; 0< T (-k 1 ) S 6 (k 2 )a^(k 3 ) - h.c. 
where 
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77o(ki, k 2 , k 3 ; I) = A)Oi, w 2 , w 3 )0 O (ki, k 2 , k 3 ; /) , 

^i(ki,k 2 ,k 3 ;i) = Di(o;i,a;2,^3)fl'i(ki,k2,k 3 ;i) , (3.20) 

with energy terms 

D (ui, u 2 , uj 3 ) = -(ui + u 2 + u 3 ) , Di(o;i, u 2 , uj 3 ) = + u 2 + io z ) . (3.21) 

The solutions of the flow equations, Eq.(3.15), for the effective coupling constants are 

g (l) = g(0)exp(-D 2 l) , 9l (l) = ff (0)exp(-D?0 , (3.22) 



which, as anticipated, eliminates the triple-gluon vertex, Eq. ( |3.17| ), for I — > oo. Corre 



spondingly, new operators in the particle number conserving sectors are generated, found 
from the second order 0(g 2 ) flow equation 

dJ! p± = iv m (l), (3-23) 

which contribute to the effective two-body, self-energy and condensate operators in the 
effective Hamiltonian. [] Solving the flow equations, Eq. ( 3.15) and Eq. ( [3.23|) , the block- 



diagonal effective Hamiltonian, H e ff, renormalized to the second order, is obtained (details 
in Appendix B) 

H eff = H d (A) + H gen (A) + 5X CT (A) , (3.25) 



where Hd(A) is the particle number conserving part of the original Hamiltonian, Eq. ([3.14 ) 



and H gen (A) includes new operators generated via perturbative elimination of H r . The cut- 
off parameter, A, regulates UV divergences in loop integrals. Renormalization is achieved 
through second order by adding the mass counterterm, 8Xct{A), 



1 One can always eliminate the particle number changing terms, appearing in second order, 

by choosing the generator r/ 2 ) = [K^,Hr ]. The flow equation reads 

^ = [^<] + ^ (2 U^], (3.24) 



and, due to the second term in Eq. ( 3.24| ), decays exponentially with flow parameter. 
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&Y OT (A) = m^ T (A)Tr J dxA 2 (x) 



m^A) = -^nJ±A* , (3.26) 

7T 6 

which absorbs the leading UV divergences in H e ff as A — ► oo. Therefore we omit the cut-off 
notation in H e ff, Eq. ( |3.25| ), because, through perturbative renormalization, the cut-off 
sensitivity is weak (see below). 

We summarize the matrix elements of H e ff in the sectors of interest (up to two gluon 
states) 2 

(0|ff e //|0) = O ren + L+C (A) + 49 (A) + O gcn (A) , 

(i|ff e //|i) = K ren + n L+c (A) + n 49 (A) + n gen (A) , 

(2\H eff \2) = V L+C + V ig + V gen , (3.27) 

where |0) is a shorthand notation for the zero-gluon sector (also vacuum state), |1) the 
one-gluon sector, etc. One should distinguish between two types of H e ff terms. The first 
one arises from normal ordering the original Hamiltonian: the instantaneous interaction 
with linear plus Coulomb potentials, labeled by L + C, and the four-gluon vertex, labeled 
by 4g. This leads to the polarization operators Hl+c and n 49 in the one-body sector, and 
condensate terms Ol + c and 4s in the zero-body sector. The energy of the ground state O 
comes from normal ordering the free Hamiltonian Hf ree with respect to the vacuum |0). As 
shown below, given the phenomenological potential Vl, the energy O of the ground state 



reproduces the known result ^from sum rules for the nonperturbative gluon condensate |2T 
The second type of terms are dynamical operators, generated by flow equations and labeled 
gen. In Eq. ( |3.27| ) the renormalized condensate O ren and kinetic K ren terms are 



°™ ^ o + = i(* - 1)V / ^ + . k ) 



2 Since calculations are done through second order in the coupling constant, new terms are gener- 
ated only up to the two-body sector (no three-body or higher interactions). 
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= (iV < " 1)v / + 2 {N * ~ 1)v / w? {— * H ' (3 - 28) 

K =K + 5X CT =j j^uw (k)a, (k) + - J — ^ uj k j 

x (af(k)<(k) + ~«(k)<(-k) + h.c.)) , (3.29) 

where O and K are defined by Eqs. (Al) and (A3), respectively, and SX^t V is the mass 
counterterm, given by Eq. (|3.26|) in the zero-body sector (analogous for 5Xq^ v ). The mass 
counterterms 5Xqt cancel the leading A 2 contribution from radiative corrections to the 
vacuum and kinetic terms. The subleading logarithmic dependence remains, producing a 
slow cut-off dependence of eigenvalues (see below). In Eq. (|3.27|) the corrections to O ren and 
K ren , regulated by the exponential cut-off function, include the condensate terms (Appendix 
B) 

Ol + c(A) = l~ Q N c (Nl - 1)V / ^^^(k " q)K " (f + {kqf) e^^ 2 , 
0, 9 (A) = ^fNM - 1)V / (3 - (k) 2 ) e-^^ 2 , (3.30) 
O gcn (A) = -a s *N c {Nl - 1)V / , ^ (k ' q) --(M-H k W k -H) 2 /A 2 , 

7 (27T) b 3w k ^ q CJ k _q (Cj k + ^ q + ^k-q) 

and the polarization operators (Appendix B) 

IWA) = \n c [ ^—V L+C (k - q)(l + (fcg) 

K + ^)«f (k)<(k) + (c 2 - ^)i«(k)<(-k) + h.c/ 



x 



n 4s(A) .^J^-i- ( 3-w,V«- 

x (V (k)<(k) + i«(k)<(-k) +h.c. 

dkdq 1 G(k, q) 



n gen (A) = -a s nN c J 



_ e -V/A 2 



(27r) 6 ^ k U; q ^ k _q Wq + LU k - c 

x (af (k)<(k) + i«(k)<(-k) + h.c.)) , (3.31) 
containing a tensor structure in the generated term 

G(k, q) = 2(1 - {kqf) (k 2 + q 2 + ^^, (1 + (L?) 2 )) , (3.32) 
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and Vl+c is the momentum space representation of the linear plus Coulomb potentials (see 
Appendix A). Also, here and through out kq = k • q. The effective gluon interaction in the 
color singlet channel is 

y f a bcf a> b 'c f ^ 1 bt (q) b ; t ( _ q)< ^ . ^ (333) 

which describes the bound state of a glueball at rest. It includes the three interactions, 
Vl+c + V^g + V gen , that define the effective Hamiltonian in the two-body sector (2\H e ff\2), 
Eq. (13~27|) . Explicitly (see Appendix B) 



V- = ^ W 0) (l " (^^J • (3-34) 
where the tensor term in the generated interaction is 



Tijfjt (k, q) = (k — q) ^njTii'^y + njUfSn' + ninjdi'ji + rii>nj'8ij — UiTijiSiij — Tii'Tij8ij> 

+ ^^(1-(^) 2 )M*y) , (3-35) 

and rij are the components of the unit vector n = (k — q)/|k — q|. All other possible 
two-body interactions that change particle number can be eliminated by the flow equations 
through second order and only contribute in higher orders. In the next section we utilize this 
effective, gluodynamical Hamiltonian, with matrix elements given by Eq. ( |3.27| ), to obtain 
and solve the gluon gap and glueball bound state equations. 

IV. GAP EQUATION AND GLUEBALLS 

Now that we have the block diagonal, effective Hamiltonian, H e ff, valid through second 
order in the coupling constant, we can nonperturbatively diagonalize the PH e ffP sector 



3 We do not consider color octet states. 
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Hamiltonian. Further, because H e ff is also renormalized, equations for physical observables 
are free from the leading UV divergences. In subsection A we first investigate the vacuum 
by formulating and solving the gluon gap equation and also calculate the gluon condensate. 
Then we address the glueball spectrum in subsection B. 



A. Gap Equation 

The gap equation allows determination of a nontrivial vacuum with gluon condensates 
and propagating quasiparticles, here gluons with a dynamical mass. There are several ways 
to obtain this equation, the most common based upon a variational principle to minimize the 
vacuum (ground state) energy. In the Bogoliubov-Valatin or BCS approach, the variational 
parameter is the angle of transformation from undressed to dressed particle (quasiparticle) 
operators. In this work the variational parameter is the quasiparticle energy, u^, which 
defines a quasiparticle basis, Eq. (|3.11| ), where the effective Hamiltonian is block-diagonal. 
Accordingly, minimizing the vacuum energy of the effective Hamiltonian 

*"=0, (4.1) 

generates the gap equation for the unknown Using Eq. ( |3.30| ) for the condensate terms 
the following gap equation is obtained 

Ji = e + m|, r (A) + \N C J ^-Lvi +C (k - q) (l + (kfjf) - ^)e-'"'^ 

+a -« N < I WTs- ( 3 " «> 2 ) ^ ~ 2a '* N < I lh— ^-e- V/A2 , (4.2) 

where the mass counterterm mcr(A) is given by Eq. ( |3.26|) and G(k, q) by Eq. ( |3.32|) . Note 



that a momentum regulating function has been introduced in the first two integrals, while 
the third is naturally regulated from the flow equations. 

The gap equation can also be obtained by using a renormalization group condition, 
namely the invariance of Hamiltonian form under the renormalization group transformation. 
Since flow equations are in this class of transformations, the renormalized effective Hamilto- 
nian should not contain off-diagonal one-body terms of the type a(k)a(— k) or aT(k)(r (— k). 
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Imposing this condition on the matrix elements of the effective Hamiltonian in the one-gluon 
sector (l\H e ff\l), Eq. (3.27) and Eq. ( |3.31| ), yields the same gap equation as above. 

There are two types of UV divergences in the gap equation: the leading, A 2 , and the 
subleading, In A, corresponding to relevant and marginal operators in the context of the 
renormalization group, respectively. Terms from the canonical Hamiltonian Hp? generate 
both types, while the confining interaction ^from H^p only generates a In A dependence. 
Renormalization is acheived by finding the counterterms associated with the canonical terms, 
since only the canonical part is renormalizable. We do not renormalize the confining inter- 
action as this would only lead to "noncanonical" divergences in each order of perturbation 
theory. The mass counterterm, Eq. (|3.26|) , cancels the leading UV divergency in the gap 
equation, Eq. (fO|), while the uncanceled part leads to a gap energy with a logarithmic 
cut-off dependence, u;(k, A). 

We have numerically solved the linearized gap equation retaining only the dominant, 
instantaneous term, with Coulomb and confining potentials, which leads to a different cor- 
responding Coulomb mass counterterm 

^t(A) = "-^A 2 , (4-3) 

instead of vr? CT which is only appropriate for the full gap equation. For fixed cut-off A = 
4GeV, the energy dispersion, u(k), is displayed in Fig. flU). The free behavior, u(k) = k, is 
recovered at high energies while for low energies a constituent gluon mass, roughly .9 GeV 
is obtained. As shown in Fig. (Q), the analytic form, w(k) = k + m(k), with a dynamical 
gluon mass, m(k), given by 

k 

mCk) = m(0)exp( ) , (4.4) 

K 

reproduces the numerical solution using parameters m(0) = 0.9 GeV and k = 0.95 GeV. 
The sensitivity of the effective gluon mass at zero momentum, m(0), to cut-off A is displayed 
in Fig. (|). 

The gluon condensate is another nonperturbative quantity calculated within this ap- 
proach. The condensate is given by the vacuum expectation value O = (0|II 2 + B^|0), Eq. 
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(3.28), but not renormalized by a counterterm. We regulate this value by subtracting the 
perturbative contribution with cu(k) — k yielding 

( i W vr3 Jo s 2w(k) • 1 j 

Using the above dispersion relation, u;(k), and a cut-off A = 4 GeV that we used in our 
previous work 0, the computed gluon condensate is .013 GeV 4 in agreement with sum rule 



results |[21| . The condensate cut-off sensitivity is shown in Fig. (|J). It is interesting that 
using the same cut-off as in 0, which calculated a condensate of .012 GeV, produces a similar 
result even though the two calculations are quite different. The logarithmic dependence of 
the effective gluon mass and condensate can both be absorbed by the running coupling 
constant a s (A). 



B. Glueball Bound State Equation 

We model the glueball bound state as two valence constituent gluons. The glueball wave 
function in the rest frame is 

M = J ^^'(q)af (q)<(-q)|0> , (4.6) 

where repeated indices are summed and a^(q) creates a quasiparticle from the nontrivial 
vacuum |0) with energy u)(q). Since in this constituent basis our effective Hamiltonian H e ff 
is block-diagonal, mixing with three gluon and higher states is suppressed. Therefore we 
can directly diagonalize in the two gluon P space and our bound state equation follows by 
projecting the fundamental equation H e ff\ip n ) = E n \ip n ) on the P space 

<Va# e// ,a?(q)af (-q)]|0) = (E n - £o)^,(q)0?'(q) , (4.7) 

where 

x^(q)0?'(q) = Ai'(q)%'(q)0?"(q) + A^q)^v(q)#?'(-q) • (4.8) 
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The polarization term, Dij, is given by Eq. (|3.13|) and H e ff by Eq. ( |3.25| ), with matrix 



elements given in Eq. Q3.27I ). Equation ( [4.7|) determines the glueball mass, M n = E n — E , 
where E is the vacuum energy defined by i? e //|0) = E \0). 

Using our H e ff we have evaluated Eqs. (4.7,4.8) and summarize in Appendix C the 
complete form of the bound state equation with all possible terms through second order. 
As mentioned in the introduction and important to repeat, even though our analysis has 
the appearance of the TDA, we do not make this approximation (truncation to the two 
gluon sector). Quite the contrary, we include the truncation corrections which introduces 
complicated new terms in the bound state equation as detailed in Appendix C and below. 

As a numerical application we consider here only the dominant, instantaneous part of the 
kernel, neglecting perturbative transverse gluon exchange and terms ^from the four-gluon 
vertex. The bound state equation for glueball states having total angular momentum, J, 
parity, P, and charge conjugation, C, is then 

" 'q 2 + m 2 CT {A) 



M n (j) n (q) 



0n(g) 



- Ji N ° / T^s WP - q) K + " q) V g (p, OUP) ■ (4.9) 

O J {271)° UJpUq 

The spectroscopic terms, F JPC , for the scalar, ++ , and pseudoscalar, 0~ + , states are 

F°++(p,q) = l + (M) 2 , 

F°- + (p,q) = 2(pg). (4.10) 



Here the Coulomb counterterm, mcr(A), is given by Eq. ( [4.3| ) and Vl+c is the sum of linear 
and Coulomb potentials in momentum space. 

Besides the UV divergences from the Coulomb interaction, the bound state equation, 
Eq. fl4.9Q , also has an IR singularity from the confining potential. The UV divergences in 
the first integral are regulated by the exponential cut-off function and the leading divergent 
part is canceled by the mass counterterm rhcri^)- The UV behavior of the second integral 
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is regulated by the wave function, which decreases with increasing momenta. Although the 
kinetic (self-energy term) and potential parts both contribute IR divergent pieces, the bound 
state equation is IR finite due to complete cancelation for small momenta. This cancelation 
only occurs for the color singlet state | fL9|j . 

Numerical solutions of Eq. ( |4.9| ) are obtained variationally with a set of gaussian test 
functions. Results for the ground state scalar and pseudoscalar glueball masses, and first 
excited states, are presented in Table | and compared with recent quenched lattice data 



22|, generated with an improved SII action and anisotropic lattice. It is more consistent to 



compare with quenched lattice results, rather than improved unquenched treatments, since 
we also omit the quark sector. Notice agreement is better in the scalar channel and that the 
mass of this state is about twice the constituent gluon mass. 



V. CONCLUSIONS AND OUTLOOK 

In this work we have described in detail a systematic procedure for solving a bound state 
problem in QCD. Our study is motivated by the success of the constituent quark model. 
Central to our approach are the flow equations, which we have applied perturbatively to 
a QCD motivated Hamiltonian with confining interaction. The flow equations generated a 
block diagonal, effective Hamiltonian which provides an improved framework for vacuum and 
excited hadron state investigations. In the block diagonal process, the effective Hamiltonian 
was also renormalized through second order in the coupling constant. In particular, we 
determined the appropriate mass counterterm to cancel the leading UV divergences. 

Applications to the gluon vacuum and excited glueball spectrum produced a nonlinear 
gap equation and a bound state equation that is superior to the TDA. Numerical solu- 
tions reproduced the QCD sum rule gluon condensate value and also provided reasonable 
agreement with quenched lattice glueball calculations. Our results support the constituent 
picture for gluonia. 

While we have only considered pure gluodynamics in this work, the same approach can 
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be directly applied to the full QCD Hamiltonian, including dynamical quarks, where the 
issue of chiral symmetry breaking is present. Of keen interest will be the structure of the 
full gap and bound state equations in the combined quark and gluon sectors. In particular, 
application to exotic four quark and hybrid systems will provide significant new information, 
especially the degree of mixing between qq, qqqq, qqg and gg states. Much of this work is in 
progress and will be reported in a future publication. 
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APPENDIX A: COMPLETE QCD MOTIVATED HAMILTONIAN 



Here we specify, in second quantized form, the full QCD motivated Hamiltonian of glu- 
odynamics in the Coulomb gauge, decomposed in the basis Eq. ( j3.11| ). One-body operators 
and condensate terms arise from normal ordering with respect to the trial vacuum state |0). 
The upper index over a Hamiltonian operator (e.g. H^) denotes the order (power) 

in (bare) coupling constant. 

Free gluon part: Eq. ( ^.4[ ) includes the gluon kinetic energy 
dk 



2 



(2k) 



;- + ^K f (kK(k) + _ c k )I (< (k)<(-k) + h. c .) 

u; k u; k 2 



/^ k af(k)<( k) + i/^(^-. k ) 



(27 



x (of (k)<(k) + i«(k)<(-k) + h.c.)) = + (K- AW) 



with 



K (0) 



dk 



■u^af (k)<(k) 



(Al) 



(A2) 
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and condensate term 



dk ,k 2 



dk k 2 



(2vr) 3 V 



( cu k ) = O (0) + (O - O (0) ) , (A3) 



with 



(A4) 



and volume V = (2tv) 3 5^(0). 

Instantaneous interaction: Eq. (|3.10 ) includes the linear confining and Coulomb 
interactions 



cabc rade 

L+C— Q J J 



x: 



rr(2) 



a'CkO+aJ'C-kO 



n^)w 3 ^(EM(^) 1/2 wk 1+ k 2 



affk-O -aff-k 



(A5) 



aj(k 3 ) + of (-kg)J ^ e (k 4 ) - af (-k 4 ) 
where the supercript (2) refers to the Coulomb term, which is second order, and Vl+cQ<.) is 

(A6) 



V L+C (k) = 2nC ad ~ + An 



Terms arising from normal ordering are the one-body operator 
N c f dk dq 



n L+ c = — 



A J (2nf 
N c f dk dq 
j2nY 



V L+C (k + q) 
V L+C (k + q) 





Da 








A j 


(q) 


V ^k^q J 



(A7) 



and the condensate term 



Ol + C = - 1)V / ^V L+C {k + q) ( ^ - 1 ) (1 + (kqf) . 



(2tt)< 



(A8) 



Nonabelian gluon part: Eq. ( |3.7| ) includes to order 0(g) the triple-gluon coupling 

(A9) 



H (l) = _W_ jabc 



2^/2" 



x: 



'n^) (2^(^k m )-^= 

Kn=l ( 27r ) / m y/v&Ws 



?(k0 + afC-kOl [aj(k 2 ) + af(-k 2 )l [<(k 3 ) + af (-k 
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where u% = u^, etc. In order 0(g 2 ), the normal ordered four-gluon vertex is 



H 



(2)_ Q-s7f rabc fade 



r bc r de f ( n 

J \n=l 



dk r , 



1 



x: 



49 A J J J Vn=l ( 27r )V ' 7 ^ "'" / v^I^iZJI 

4(k0 + af (-kO] [ S c (k 2 ) + af (-k 2 )] U d (k 3 ) + af (-kg)] k e (k 4 ) + af (-k 4 



Also, after normal ordering one obtains a one-body operator 

dkd<\ 1 



„ Ar f dkdq 1 > , 1 

n 4, = «^ Cj / [2^-A,(q)] 



X 



(27r) 6 WkCUq 

<t(k)a J a (k) + i(<(k)a»(-k)+h.c. 



and a condensate term 



o 49 = ^a c( a c 2 -i)v ^ krfq 1 



3 - (kqf 



4 " c ■ y (2vr)6^ q 
One-body and condensate terms diverge in the UV region and must be regulated. 



(A10) 



(All) 



(A12) 



APPENDIX B: SECOND ORDER FLOW EQUATIONS 

Solving the second order flow equation for the particle number conserving part Hd, Eq. 
( p.23| ), generates three types of Hamiltonian operators: two-body effective interactions, one- 
body polarization terms and condensates. We consider each separately. 



1. Effective Interaction (two-body sector) 

We calculate an effective gluon interaction in the color singlet channel for glueball Fock 
states. The general form of an effective interaction is given by Eq. ( |3.33| ) and includes 
interactions from the original Hamiltonian and a new interaction, generated to second order 
by flow equations in the two-body sector. Consider first the generated terms. The t-channel 
diagrams, arising from the commutator [a^a^a^., ai4' a 2V a 3'fc / ] (diagram without backward 
motion) and from the commutator [aij£4j a 3fcj a i'l' a 2'j' a 3'l'] (Z-graph in t-channel), do not 
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contribute to the color singlet state; only the s-channel diagram, coming from the commu- 
tator [ctiJa2j a 3A,o a i'j' a 2'j' a 3'l']' nee ds to be calculated (here the notation is = a®(ki), etc.). 
In the cm. frame the flow equation for the effective interaction is 

dV gen (l) I £abc£a > b > c fdkdq 1 



j?abc ja'b'c J 



dl 8 J (2ir) 6 CUkCUqCUk-q 

x r iifc (k, -q, -(k - q))r i y A .,(k, -q, -(k - q))D kk ,(k - q) ■ 4 • 9 

x foi(q, k, k - q; l)g x (k, q, k - q; Z) + »7i(k, q, k - q; %i(q, k, k - q; /)] 

x :af(q) a f(-q)<(kK'(-k):, (Bl) 



where the factor 4 is the number of permutations and the property of T-factors, Eq. ( |3.18[) 



is used. The generators are given in Eq. ( p.20[) and coupling constants in Eq. (|3.22| ). The rjg 
sum corresponds to two different time-ordered s-channel diagrams. We introduce the factor 
Sijk,i'j'k'(K q) by 

rij- fc (k, -q, -(k - q))r i / J -/ fc /(k, -q, -(k - q)) = ~4S , ii fc ii # J v Jfe /(k, q) . (B2) 

which, due to transversality, can be expressed as 

Sijk,i'j'k'{K q) = (kj5 ik - ]-{k + q) k 5 i:j + qi8j k ){kj'8 V y ~\{k + q)k>$i>j> + qi>8j>k>) , (B3) 

and is symmetric under interchange of k and q. The tensor structure of the generated 
interaction is defined by the contraction Sijk,i>j'k>(k, q)Dfcfc'(k — q) = T^y (k, q), where 



Tij^jiik, q) = [qiqi>5jji + kjkj>8u> + qi{kj>8 V j - kj8i>j>) + qi>(kj5ij> - kfSij) 



(k — q) 2 iuiUi'Sjjt + njUj>5ii' + UiUjSi/j' + ThitUj/Sij — uinybiij — ni>rij5i 



and are the components of the unit vector n = (k — q)/|k — q|. Integrating the flow 
equation, Eq. (|HT|) , we obtain the generated interaction V gen {l — > oo) = V gen (the initial 



value is V gen (l = 0) = 0) 
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dkdq 1 _ .D 1 + D v 



V gen = a s 2^r bc f a ' b ' c I =f = T^ Y (k,q) 



x :af(q}ap(- q )a^k)^(-k):, (B5) 

with 

Dx = -(-LO k + ujq + w k _ q ) , L>r = -(-w q + + ^k-q) • (B6) 
This reduces to 



„!„■ ,„'//,■ I dkdq 1 / U k - w q ) 



2 



^ gen = -a s 27rr 6c r 6c / ttttI 2" ?W(k,q) 1 



x ^J^^a^VqXCkX^-k):. (B7) 

Now consider the second order gluon interactions that conserve particle number. The 
generator Eq. ( 3.19|) , produces second order interactions that change the number 



of quasiparticles, for example, an interaction W having form a^a^a^a^. However, the 
generator choice, r] = f/^+r/ 2 ), eliminates W, where the second order generator is calculated 
from r/ 2 ) = [K^°\W}. 

To maintain gauge invariance we also include the normal ordered instantaneous, Eq. 
(K5|), and the four-gluon, Eq. ( |A10|) , interactions of the original Hamiltonian in the color- 
singlet channel. Projecting these terms on the color singlet state yields 

V L+ c = -\f*°f«"f ^ ( ^ q) V L+c (k - q) M*V :af(q)a;; + (-q)<(k)<;(-k): 
V i9 = ^f^f^J (6 ul 8 jf - :af(q)ap(-ci)at(k)4(-ky. , (B8) 

The complete effective gluon interaction is given by 

V eff = V L+c + V 4g + V sen , (B9) 
with the generated interaction defined by Eq. 



2. Polarization Operator (one-body sector) 

Two contraction terms from the commutator fr/ 1 -*, H^J] contribute to one-body operators 
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n = / ^3 (n« 6 (k)af (k)aj(k) + Mj*(k) i (a?(k)a$(-k) + h.c.)) , (BIO) 
with polarization amplitudes, nfj'(k), (k), given by gluon loop integrals of the form 

n ^ (k)= /(^^ (k ' q) - (B11) 

The specific integrand structure is given below. Note, from rotational invariance, the one- 
body operators n°*(k) and M^ 6 (k) have tensor structure, A ab (k)5ij + B ab (k)kikj, however, 
due to transversality, only the term proportional to Sij survives after integrating Eq. ( |B10| ) . 
Hence, we only need to extract A ab (k) which can be achieved by choosing k = k z and 
considering components i = j = x and y but 7^ z. Then in Eq. ( |B11| ), where 11°* (k, q) also 
has terms of the type 5^ and qi<]j (there are no terms proportional to hi and kj), one can 
substitute 

m^\{i-{k) 2 )Sa, (B12) 



to extract the 5ij component in Eq. (|B 1 1| ) . ^From the second order flow equation, Eq. 



( [3.231 ), the gluon polarization operator is 



dU(l) 1 ^ f dkdq I 



dl 8 J (27r) 6 dJkdJqtUk-q 

x Y ijk {k, -q, -(k - q))r i y A .,(k, -q, -(k - q))D jjl (q)D w (k - q) • 2 • 9 • 2 

(771 (k, q, k - q; %l(k, q, k - q; /) + ?7o(k, q, k - q; l)g (k, q, k - q; I)) af (k)«#'(k) 
+ (Vl (k, q, k - q; l)g (k, q, k - q; /) + ?7 (k, q, k - q; %l(k, q, k - q; I)) (B13) 
x i«(k)<;(-k) + h.c.) , 

where r\\g\ corresponds to a gluon loop without backward motion, 770^0, a Z-graph, and the 
last two terms rjigo and rjogi, are gluon loop diagrams with two incoming and two outgoing 
gluon lines. In calculating Eq. ( |B13| ) we have also used the property, Eq. ( |3.18j ), for 
T-factors. 

The polarization operator tensor structure is given by a double contraction of the S^k, q) 
factor, Eq. ( |B3|) , with polarization terms in Eq. flB13|) . We then define 
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F u >(k, q) = .S, ;/.,,' ; v/fk. q}Dj jl (q)D kk/ (k - q 
11 
<1 



5u>k 2 q 2 ( {- 2 + /u - {kqf) ) + 2g^ ( 1 



A: 2 



r (l - {kqf) , (B14) 



(k-q)^- — 7 ' — ^ 2(k-q)»' 

with repeated indices summed. Due to transversality, only the 8^/ component is retained, 
therefore from Eq. (II 



F ft ,(k,q) -><5«,G(k,q)/2, 
G(k,q) = 2(l-(^) 2 ) (V- <? + 



k 2 n 2 

2,1,2 2 Ay (1 + M 2j 



(B15) 



2(k-q) 2 

We integrate the flow equation, Eq. (|B13|) , with the generators and coupling constants given 
by Eq. ( p.20|) and Eq. (|3.22|) , respectively, producing the second order polarization operator 
correction, 511 = 11(7) — H(Iq = 0). The polarization operator, at a scale I = 1/A 2 , is then 

dkdq 1 



n g en(A) = a s irN c J 



(27r) 6 UkUJqUJk-, 



-G(k,q) 



X 



+ 



/ 1 



\2D 



2D X 



a- T (k)<(k) 



^^e-(-o^)A^ K(k ) <( _ k ) +h . c .) 



(B16) 



with 



Dr 



(B17) 



For large momenta flowing in the loop, one has D ~ D\ ~ — (u; q + co>k_ q ) ~ — 2cj q . Finally, 
rescaling the cut-off, A — > \/2\, which does not change results, the polarization operator 
generated by the flow equations through second order is 

dkdq G(k, q) 1 



n gen (A) = -a s 7iN c J 



(27r) 6 w k u; q u; k _ q ^ q + 



-4g 2 /A 2 



x (af(k)<(k) + ^«(k)<(-k)+h.c.)) , 



(B18) 



where in the exponential factor we have used the free dispersion relation, u; q = q, valid 
for large cut-off values. We then calculate one-body operators, by normal ordering the 
instantaneous interactions and the four-gluon vertex, Eq. (A5) and Eq. ( A10|) , respectively, 
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and again extract the 5ij component using Eq. (P12| ). The regulated polarization operators 
are 

IWA) = \n c I ^—V L+C (k - q)(l + fa)*)**'* 

8 J {271)° UJkUJq 

x (K + u£)of (k)<(k) + ( W J - ^)l«(k)<(-k) + h.c.)) (B19) 
x faf(k)<(k) + l«(k)<(-k)+h.c.; 



2 

We choose the same regulating function to match energy denominators. Note that in the 
generated term the regulator emerges naturally /,from the flow equations, while the normal 
ordered terms require introducing a loop integral regulator with form consistent with the 
generated term. The complete polarization operator is then 

n(A) = iw(A) + n 4 ,(A) + n gen (A) . (B20) 

We conclude this subsection by specifying the counterterm for the leading UV divergences 
in the polarization operator. The instantaneous (Coulomb only), four-gluon and generated 
terms all contribute to the quadratically UV divergent part, H dw , with respective weights, 
1,2,-1/4 

n^(A) = ^A 2 (l + 2-l) 
7r 3 V 4/ 

x / 7^1^ (< f (k)<(k) + ~«(k)<(-k) + h.c.)) . (B21) 
The appropriate counterterm is then 

dk m 2 CT (A) (a\ n ,a n , , 1 



(2vr) 3 2u; k 
m 2 CT (A) 



5X' CT {A) = J (2n)3 mC 2 l[ } (af(k)<(k) + -«(k)<(-k) + h.c.) 

'dx:A°(x)A a (x):, (B22) 



with mass 

m 2 CT (A) = ~-N c ^A 2 . (B23) 

7T O 

When the quark sector is included using this same procedure, the counterterm coefficient 
reproduces the QCD /5-function. This result was also obtained in our previous study ||. 

31 



3. Gluon Condensate (zero-body sector) 



The second order flow equation, Eq. ( |3.23j ), for the condensate term is 
dO gen (l) _ l ATfAT 2 1V ,, f dkdq 1 



-N C (N 2 C - 1)V J ■ 



dl 8 C 7 (27r) 6 CJkCJqto'k-q 

xr^ fc (k, -q, -(k - q))r i / 3 -/ fe /(k, -q, -(k - q))D i i'(k)I> i ,y(q)I? fcfc /(k - q) • 6 • 2 
X7? (k, q, k - q; l)g (k, q, k - q; I) , (B24) 

where only the commutator [a>xi a 2j a 3k> a i4' a 2'i ,a 3'l'] contributes to the vacuum expectation 
value. The factor of 6 is the number of permutations and the volume is V = (27r) 3 5^ 3 ^(0). 

The tensor structure of the condensate term is given by the triple contraction of the 
5-factor, Eq. (p3|), with polarization terms, i.e. Sijk^j'k'iM, <l)Diii(k)Djji(q)Dkki(\i _ q) = 
G(k, q), where the function G(k, q) is given by Eq. ( |B15|) . Integrating yields the correction, 



50 = 0(1) — O(lo = 0), where 0(1) is a condensate for flow parameter I, related to the 
energy scale A by I = 1/A 2 . The resulting generated condensate term through second order 
is 

0,„(A) = a,*N c (N? - 1)V / ggjS 1 y ,-^ , (B25) 

J (2n) b w k a; q u; k _ q 3D 

where the energy difference Dq is given by Eq. ( |B17| ). ^ The three regulated condensate 



terms, instantaneous, Eq. (|A8|) , four-gluon, Eq. (|A12|) , and generated can be summarized 



4 In the gap equation calculations we use the symmetrized form 

dkidk2<ik3 1 



Ogen(A) = a s ltN c (N 2 _ 1) V j 



(2tt) 9 UJ1LU2UJ3 



x (2vr) 3 ^ 3 ) ( kl + k 2 + k 3 ) G(kl * 2 > k3) e- 2 ^/A 2 

SDq 



where -Do is given in Eq. ( 3.21 ), and 

G(k!,k 2 ,k 3 ) = [(k 2 + k 2 + k 2 )(i-2(kk)(kh)(kk)) 
+k 2 (k 1 k 2 ) 2 (k 1 h) 2 + *${kk) 2 {kh) 2 + *$(kh) 2 (hh) 2 ) 

X (2vr) 3 <5( 3 ) (ki + k 2 + k 3 ) G(k, q) . 
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O l+c {\) = \NM - 1)V / ^V(k - q) - l) (1 + (k) 2 ) e-'^' 

a ' n ATfAT2 i^7-/" dkd< l 1 ^ /-£^2> „-(|g|+[fc|) a /A a 



4g (A) = -^iV c (iV c 2 - l)V / -— J 3 - e-MW/A- (B26) 



O gcn (A) = -a s vriV c (iV c 2 - 1)V | 



dkdq 1 G(k, q) c -(M+|fc|+|fc-g|) 2 /A 2 

(27r) 6 3^ k ^q^k-q (^ k + ^ q + ^ k -q) 



The regulating procedure is the same as above. The complete radiative correction to the 
gluon condensate, Eq. (|A3|), is then 



0(A) = L+C {\) + 49 (A) + O gcn (A) . (B27) 

Finally, we determine the counterterm for the leading UV divergences in the zero-body 
sector. Again, the instantaneous (only Coulomb part), four-gluon and generated terms 
contribute, respectively 

ldiv/AN « s 2iV c (iV c 2 -l)V A2 ^ io 1\ f dk 1 



O dlv (A) = — — ^ — A 1 + 2 - - / -—- . (B28) 

V ; 7T 3 V A) J (27r) 3 2c;k V ' 



Hence the counterterm is 



8X>> T (k) = {Nl - 1)V / ^ ! %^ = ! ^ / rfx(0|A a (x)A a (x)|0) , (B29) 
where the mass mcr(A) is given by Eq. ( B23|) . Combining the counterterms from one-body, 



Eq. ([B22|) , and zero-body, Eq. (|B29| ), sectors yields 

SX CT = 5X' CT + 5X'c T = m 2 7T (A)Tr J M 2 (x) , (B30) 

which is the complete mass counterterm in the effective Hamiltonian, renormalized through 
second order, Eq. ( |3.25| ). Note, that the one- and zero-body counterterms have the same 
mass coefficient, mcr(A). This suggests one can work directly in a field theoretical basis 
without a Fock operator decomposition. 
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APPENDIX C: COMPLETE GLUEBALL BOUND STATE EQUATION 



In this appendix we present the complete bound state equation, including all second 
order terms. Projecting the effective Hamiltonian, Eq. fl3.25| ), on the two-gluon Fock basis, 
Eq. ( |4.7| ), generates the bound state equation for a glueball at rest with mass M = E — Eq 
(we suppress the excited state quantum index number) 

'q 2 + m 2 CT {A) 



M(^(q)0f(q) + ^'(q)^(q) 



dp 



+ \ N 'J (ar). 



+ a s 7riV, 



V L+C (p - q) (1 + 
dp — (3 - Cm 2 ) e" p2/A 



■2\ U P + U; q c -p 2 /A 2 



(27r) 3 CUp^ q 

dp 1 



G(p, q) 



^-4p 2 /A 2 



(27r) 3 w p u; q c i ; p _ q cj p + ^ p _ q 

(jfS(q)^(q) + i£(q)^(q); 



+ 

+ a s 7riV, 



+ o> s 2ttN c J 



(2tt) 3 
dp 2 
(27r) 3 a;pa; q 
c/p 1 



iv B / ^ Wp - q) ^ . ' , 7 M/ (^L(q)^n m (p)^ w (p) + n^(q)^ m (p)0'- n (p) 



Ai(q)An(p)^(p) - ^(q) [Xf„ m (p)C(p)+C(p)C(p) 



(^ p -^ q ) 2 + u;p_ q/ 



x 



(2?r) 3 WpW q W^ q \ 

T fcm , fc w(q,p) (xg f (q)xir'(p)^(p) +^(q)^r'(p)^ n (p: 



(CI) 



where the factor T^ m ^i m i is defined in Eq. (B4) and 



>£(p) 



(.in , 



p)±<T(- P )) 



Xg(q) = Ai(q)^n(q) + An(q)^((q) 
K^'(q) = £>«(q)D in (q) - An(q)%(q) 



(C2) 



involving the polarization sum Z)j 3 (q), Eq. (3.13), and glueball momentum wavefunction, 
(p ln (p), introduced in Eq. (4.6). 

The bound state equation can be further reduced by forming the appropriate tensor 
contractions for u (p) corresponding to a specific glueball state having quantum numbers 
J, total angular momentum, P, parity, and C, charge conjugation. For the applications in 
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this paper we considered the scalar (J PC = ++ ) and pseudoscalar (J PC = CT + ) channels 
having tensor contractions y '(q) = S' lj (f)(q) and lJ (q) = e tjk q k (j)(q), respectively. The 
resulting equation describing both J = states is 

" 'q 2 + m 2 CT {A) 



Mc/>(q) = 



2\ ^P + ^q c -pVA 2 



+ a s irN c 



dp 1 

(27r) 3 CJpCJq 

dp 1 



(3 - (kq) 2 ) e~ p: 



G(p, q) 



(27r) 3 a; p a; q a;p-„ q u; p + uj p - 



/a 2 



,-V/a 2 



+ 



4* 



+ a s irN c 



dp 

(2^)3 
' dp 



(2tt) 3 2a; p a; ( 



Vi+c(p - q) 



"(3 



[Up + CJ q ) jpc- 



2 )£ JPC (p, q)0(p) 



F^(p,q)0(p) 



+ a 



7 (27T)' i CJ p U; q Wp_ q \ 



(w p - u n f 



(iOp - CJ q ) 2 + 



p-q, 



with spectroscopic terms D JPC , E JPC and F JPC for each channel given by 



(C3) 



F 0++ (p,q) = l + (p^ 2 
F°- + (p,q) = 2(M) 
£°++(p,q) = l 
E°- + (p,q) = 



D u++ (p, q) = G(p, q) = 2(1 - (pqY) p 2 + g 2 + 



D u - + (p,q) = 2pV(l-(j3g) 2 ) - + 



p 2 q 2 



2(p-q) 5 



(C4) 



\P? (p - q) 2 

An approximate form for Eq. (C3) is used (see Eq. (4.9)) in the main text for numerical 
calculations. 
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TABLES 





0++ 


0*++ 


o-+ 


*-+ 


calculated (MeV) 


1760 


2697 


2142 


2895 


lattice data H (MeV) 


1730(80) 


2670(130) 


2590(130) 


3640(180) 



TABLE I. Glueball spectrum for the first two scalar and pseudoscalar states 
(a 5 = £ = 0.4, a = 0.18 GeF 2 , A = AGeV). 



38 



FIGURES 




1 2 3 4 5 6 

q(GeV) 

FIG. 1. One particle dispersion relation, w(q). Dots are the numerical solution of the gap 
equation for parameters a s = 0.4, a = 0.18 GeV 2 and A = AGeV . The solid line is the free 
dispersion, w(q) = q. 



m(GeV) 




1 2 3 4 5 6 

q(GeV) 

FIG. 2. Dynamical gluon mass. Dots represent the numerical solution for m(q) = cj(q) — q 
(same parameters as in Fig. 1). Solid line is m(q) = 0.9exp(— g/0.95) (parameters are in GeV). 
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FIG. 3. Cut-off dependence of the constituent gluon mass (same parameters as in Fig. 1). 
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FIG. 4. Gluon condensate cut-off dependence (same parameters as in Fig. 1). 
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